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Abstract  The  present  paper  highlights  results  derived 
from  the  application  of  a  high-fidelity  simulation  technique 
to  the  analysis  of  low-Reynolds-number  transitional  flows 
over  moving  and  flexible  canonical  configurations  moti¬ 
vated  by  small  natural  and  man-made  flyers.  This  effort 
addresses  three  separate  fluid  dynamic  phenomena  relevant 
to  small  fliers,  including:  laminar  separation  and  transition 
over  a  stationary  airfoil,  transition  effects  on  the  dynamic 
stall  vortex  generated  by  a  plunging  airfoil,  and  the  effect  of 
flexibility  on  the  flow  structure  above  a  membrane  airfoil. 
The  specific  cases  were  also  selected  to  permit  comparison 
with  available  experimental  measurements.  First,  the  pro¬ 
cess  of  transition  on  a  stationary  SD7003  airfoil  section  over 
a  range  of  Reynolds  numbers  and  angles  of  attack  is  con¬ 
sidered.  Prior  to  stall,  the  flow  exhibits  a  separated  shear 
layer  which  rolls  up  into  span  wise  vortices.  These  vortices 
subsequently  undergo  spanwise  instabilities,  and  ultimately 
breakdown  into  fine-scale  turbulent  structures  as  the 
boundary  layer  reattaches  to  the  airfoil  surface.  In  a  time- 
averaged  sense,  the  flow  displays  a  closed  laminar  separa¬ 
tion  bubble  which  moves  upstream  and  contracts  in  size 
with  increasing  angle  of  attack  for  a  fixed  Reynolds  number. 
For  a  fixed  angle  of  attack,  as  the  Reynolds  number 
decreases,  the  laminar  separation  bubble  grows  in  vertical 
extent  producing  a  significant  increase  in  drag.  For  the 
lowest  Reynolds  number  considered  ( Rec  =  104),  transition 
does  not  occur  over  the  airfoil  at  moderate  angles  of  attack 
prior  to  stall.  Next,  the  impact  of  a  prescribed  high-fre¬ 
quency  small-amplitude  plunging  motion  on  the  transitional 
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flow  over  the  SD7003  airfoil  is  investigated.  The  motion- 
induced  high  angle  of  attack  results  in  unsteady  separation 
in  the  leading  edge  and  in  the  formation  of  dynamic-stall- 
like  vortices  which  convect  downstream  close  to  the  airfoil. 
At  the  lowest  value  of  Reynolds  number  ( Rec  =  104), 
transition  effects  are  observed  to  be  minor  and  the  dynamic 
stall  vortex  system  remains  fairly  coherent.  For  Rec  =  4  x 
104,  the  dynamic-stall  vortex  system  is  laminar  at  is 
inception,  however  shortly  afterwards,  it  experiences  an 
abrupt  breakdown  associated  with  the  onset  of  spanwise 
instability  effects.  The  computed  phased-averaged  struc¬ 
tures  for  both  values  of  Reynolds  number  are  found  to  be  in 
good  agreement  with  the  experimental  data.  Finally,  the 
effect  of  structural  compliance  on  the  unsteady  flow  past  a 
membrane  airfoil  is  investigated.  The  membrane  deforma¬ 
tion  results  in  mean  camber  and  large  fluctuations  which 
improve  aerodynamic  performance.  Larger  values  of  lift 
and  a  delay  in  stall  are  achieved  relative  to  a  rigid  airfoil 
configuration.  For  Rec  =  4.85  x  104,  it  is  shown  that  correct 
prediction  of  the  transitional  process  is  critical  to  capturing 
the  proper  membrane  structural  response. 


1  Introduction 

Unsteady  low-Reynolds-number  flows  are  found  in  natural 
flyers,  as  well  as  in  small  unmanned  air  vehicles  and  micro 
air  vehicles  (or  MAV’s)  due  to  the  relatively  small  size  and 
low  flying  speeds  involved  (Mueller  1985;  Shyy  2008; 
Ellington  et  al.  1996).  Depending  upon  the  specific  condi¬ 
tions,  these  flow  fields  may  be  characterized  by  extensive 
regions  of  laminar  flow,  by  the  onset  of  laminar  separation 
bubbles  (even  at  moderate  incidence),  and  by  laminar- 
turbulent  transition  zones.  For  the  case  of  flapping  wings,  as 
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well  as  for  severe  gusts,  the  highly  unsteady  forcing  induces 
the  formation  of  dynamic- stall  and  leading-edge  vortices 
whose  evolution  and  interaction  with  the  aerodynamic  sur¬ 
faces  have  a  significant  impact  on  flight  stability  and 
performance.  Although  much  has  been  studied  and  learned 
about  these  unsteady  vortical  flow  features,  challenges  still 
remain  in  understanding  fully  their  structure,  scaling  and 
implications  on  flight  efficiency,  in  particular  over  the  broad 
range  of  parameters  encountered. 

Vehicle  weight  considerations  in  MAV  design  dictate 
the  use  of  lightweight  and  highly  flexible  or  compliant 
structures.  These  characteristics  create  additional  com¬ 
plexities  due  to  the  rich  fluid- structure  interactions 
generated  by  coupling  of  aerodynamic,  inertial  and  elastic 
forces.  There  is  growing  recognition  that  wing  deforma¬ 
tion  may  be  beneficial  to  flight  performance  including 
stall  delay  (Smith  and  Shyy  1995;  Shyy  and  Smith  1997; 
Song  and  Breuer  2007;  Rojratsirikul  et  al.  2008;  Waszak 
et  al.  2001).  However,  harnessing  these  passive  effects 
through  aeroelastic  tailoring  clearly  demands  further 
investigation. 

From  the  perspective  of  analysis  and  simulation,  this 
non-traditional  low-Reynolds-number  aerodynamic  regime 
over  flexible  or  flapping  surfaces  poses  a  severe  challenge 
due  to  the  following  key  factors.  Difficulties  arise  due  to 
the  presence  of  highly  unsteady  flows  which  defy  standard 
quasi-steady  characterization.  The  flow  fields  are  of  a 
mixed  laminar-transitional-turbulent  type  which  high- 
Reynolds-number  analysis  tools  are  not  particularly  well- 
suited  to  handle.  Both  in  nature  and  in  MAV  applications, 
an  extensive  range  of  parameters  and  configurations  are 
encountered.  Lastly,  there  exists  a  strong  coupling  of  the 
unsteady  aerodynamics  and  structural  response  which 
requires  advanced  multidisciplinary  approaches. 

Given  the  aforementioned  difficulties,  a  heirarchy  of 
increasingly  more  complex  canonical  model  problems  can 
be  considered  in  order  to  facilitate  progress  in  the 
improved  understanding  and  prediction  of  the  multi-dis¬ 
ciplinary  physics  relevant  to  small  flyers.  The  simplest  of 
these  configurations  is  a  maneuvering  or  flexible  airfoil 
section  as  a  model  for  flapping  flight.  In  that  spirit,  the 
present  effort  addresses  the  application  of  a  high-fidelity 
simulation  technique  for  the  analysis  of  low-Reynolds- 
number  transitional  flows  over  stationary,  plunging  and 
flexible  airfoils.  The  methodology,  summarized  in  Sect.  2 
and  described  in  more  detail  elsewhere  (Visbal  and 
Rizzetta  2002;  Visbal  et  al.  2003),  is  based  on  a  6th-order 
accurate  implicit  large-eddy  simulation  (ILES)  procedure 
incorporating  a  low-pass  spatial  numerical  filtering  tech¬ 
nique.  The  high  order  algorithm  facilitates  the  proper 
representation  of  the  transition  process  whereas  the  low- 
pass  filter  provides  regularization  in  turbulent  flow 
regions  in  lieu  of  a  standard  sub-grid-scale  model.  This 


method  is  particularly  attractive  for  low-Reynolds-number 
applications  since  it  permits  a  seamless  treatment  of  the 
mixed  laminar,  transitional  and  turbulent  flow  features 
previously  noted. 

The  first  problem  considered  (Sect.  3)  is  the  detailed 
structure  of  the  transitional  flow  over  a  rigid  stationary 
SD7003  airfoil  section.  Several  experimental  and  compu¬ 
tational  studies  of  this  configuration  have  been  performed 
recently  (01  et  al.  2005;  Radespiel  et  al.  2006;  Yuan  et  al. 
2005;  Lian  and  Shyy  2006).  Although  this  canonical 
problem  represents  perhaps  the  simple  abstraction  of  low 
Reynolds  number  aerodynamics,  it  exhibits  a  complex  flow 
physics.  Of  particular  interest  here  is  the  process  of  lami¬ 
nar-turbulent  transition  over  the  airfoil  associated  with  the 
separated  shear  layer  which  is  present  even  at  modest 
incidence.  In  the  classical  time-averaged  representation, 
this  leads  to  the  formation  of  the  so-called  laminar  sepa¬ 
ration  bubble  (LSB)  which  eventually  bursts  leading  to 
airfoil  stall.  The  effects  of  Reynolds  number  (104  <  Rec  = 
pcU^/p  <  9  x  104)  and  angle  of  attack  (2°  <  a  <  14°)  on 
the  transition  process  are  investigated. 

Of  more  interest  to  biological  and  MAV  flight,  is  the 
aerodynamics  of  maneuvering  airfoils  as  a  tractable  model 
of  flapping  flight.  As  a  point  of  departure  in  this  direction, 
we  consider  in  Sect.  4  the  unsteady  transitional  flow  over  a 
plunging  SD7003  airfoil.  The  specific  flow  and  motion 
parameters  correspond  to  the  recent  experimental  and 
computational  study  of  McGowan  et  al.  (2008).  The  airfoil 
is  set  at  a  static  angle  of  attack  a0  =  4°,  and  the  reduced 
frequency  and  plunging  amplitude  are  k  =  nfc/U^  =  3.93 
and  hlc  =  0.05,  respectively.  These  parameters  result  in  a 
maximum  excursion  in  motion-induced  angle  of  attack  of 
21.5°  which  promotes  through  unsteady  separation  the 
generation  of  dynamic-stall-like  vortices  near  the  airfoil 
leading  edge.  The  transitional  behavior  of  these  vortices,  as 
they  propagate  along  the  airfoil,  is  investigated  for  two 
values  of  the  Reynolds  number  ( Rec  =  104  and  4  x  104). 

Finally,  in  Sect.  5,  we  study  the  impact  of  structural 
compliance  on  the  unsteady  flow  past  a  stationary  flexible 
airfoil  at  low  Reynolds  number.  The  configuration 
employed  consists  of  a  latex  sheet  attached  to  two  rigid 
mounts  located  at  the  leading-  and  trailing-edge  of  the 
membrane,  and  corresponds  to  the  experimental  arrange¬ 
ment  of  Rojratsirikul  et  al.  (2008).  The  effect  of  increasing 
angle  of  attack  on  the  fluid-structure  interaction  is  exam¬ 
ined  for  a  fixed  low  Reynolds  number  ( Rec  = 2.5  x  103). 
As  the  membrane  deforms  due  to  the  aerodynamic  load,  a 
cambered  geometry  is  achieved  which  modifies  the  stall 
characteristics  of  the  airfoil.  For  high  incidence,  large 
fluctuations  in  membrane  shape  occur  which  can  be  viewed 
as  forcing  of  the  separated  vortical  flow.  Limited  results  on 
the  effect  of  transition  on  the  flow  plast  the  flexible 
membrane  are  also  considered  for  Rec  =  4.85  x  104. 
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2  Methodology 

2.1  Governing  equations 

In  the  present  moving  or  flexible  airfoil  simulations,  the 
governing  equations  are  taken  to  be  the  full  unsteady  Na- 
vier-Stokes  equations  cast  in  strong  conservative  form  and 
incorporating  a  general  time-dependent  curvilinear  coor¬ 
dinate  transformation  (x,  y,  z,  t)  ->  (£,  rj,  £,  t)  (Vinokur 
1974;  Steger  1978)  from  physical  to  computational  space. 
In  terms  of  non-dimensional  variables,  these  equations  can 
be  written  in  vector  notation  as: 

0  fu\  dF  0G  dH  _  1 
0t  y/y  0<J  dr]  +  0£  Re 

where  U  =  {p,pu,  pv,  pw ,  pE}  denotes  the  solution  vector, 
J  m  0(£,  fy,  £,  t)/0(jc,  y,  z,  t)  is  the  transformation  Jacobian, 
F,G  and  H  are  the  in  viscid  fluxes,  and  FV,GV  and  Hv 
represent  the  viscous  terms.  The  full  form  of  these  terms 
can  be  found,  for  instance,  in  Anderson  et  al.  (1984). 

It  should  be  noted  that  the  above  governing  equations 
correspond  to  the  original  unfiltered  Navier-Stokes 
equations,  and  are  used  without  change  in  laminar, 
transitional  or  fully  turbulent  regions  of  the  flow.  Unlike 
the  standard  LES  approach,  no  additional  sub-grid  stress 
(SGS)  and  heat  flux  terms  are  appended.  Instead,  a  high- 
order  low-pass  filter  operator  (Visbal  and  Gaitonde  1999; 
Gaitonde  et  al.  1999)  is  applied  to  the  conserved 
dependent  variables  during  the  solution  of  the  standard 
Navier-Stokes  equations.  This  highly-discriminating  filter 
selectively  damps  only  the  evolving  poorly  resolved 
high-frequency  content  of  the  solution  (Visbal  and 
Rizzetta  2002;  Visbal  et  al.  2003).  This  filtering  regu¬ 
larization  procedure  provides  an  attractive  alternative  to 
the  use  of  standard  sub-grid- scale  (SGS)  models,  and  has 
been  found  to  yield  suitable  results  for  several  turbulent 
flows  on  LES  level  grids.  A  re-interpretation  of  this 
implicit  LES  (ILES)  approach  in  the  context  of  an 
approximate  deconvolution  model  (Stolz  and  Adams  1999) 
has  been  provided  by  Mathew  et  al.  (2003).  The  ILES 
approach  is  very  attractive  for  the  present  applications 
involving  mixed  laminar,  transitional  and  turbulent  flow 
regions. 

The  structural  model  for  the  membrane  airfoil  is  based 
on  the  approach  of  Smith  and  Shyy  (1995)  for  a  two 
dimensional  elastic  membrane  subjected  to  a  normal  force. 
The  governing  equation  in  nondimensional  form  is 

=  A/?  (2) 

where 


T  =  Eh(50  +  5)  (3) 

5  =  L  -  1  (4) 


and 


In  these  equations  ps ,  h ,  Cd,  T,  E  and  S0  are  the  membrane 
density,  thickness,  damping,  tension,  modulus  of  elasticity 
and  the  length  increase  of  the  pretensed  membrane, 
respectively. 

2.2  Numerical  procedure 

All  simulations  are  performed  employing  the  extensively 
validated  high-order  Navier-Stokes  solver  FDL3DI  (Visbal 
and  Gaitonde  1999;  Gaitonde  and  Visbal  1998).  In  this 
code,  a  finite-difference  approach  is  employed  to  discretize 
the  governing  equations,  and  all  spatial  derivatives  are 
obtained  using  a  sixth-order  compact-differencing  scheme 
(Lele  1992). 

In  order  to  eliminate  spurious  components,  a  high-order 
low-pass  spatial  filtering  technique  (Visbal  and  Gaitonde 
1999;  Gaitonde  et  al.  1999)  is  incorporated.  The  filter  is 
applied  to  the  conserved  variables  along  each  transformed 
coordinate  direction  once  after  each  time  step  or  sub-iter¬ 
ation.  For  transitional  and  turbulent  flows,  the  previous 
high-fidelity  spatial  algorithmic  components  provide  an 
effective  implicit  LES  approach  in  lieu  of  traditional  SGS 
models,  as  demonstrated  in  Visbal  and  Rizzetta  (2002)  and 
Visbal  et  al.  (2003). 

For  the  case  of  a  maneuvering  airfoil,  the  grid  is 
moved  in  a  rigid  fashion  following  the  prescribed  kine¬ 
matics  which  provides  the  position  and  velocity  of  the 
mesh  at  every  instant  in  time.  For  the  flexible  membrane 
airfoil  computations,  the  aerodynamic  mesh  is  allowed  to 
deform  in  accordance  with  the  motion  of  membrane.  A 
simple  algebraic  method,  described  in  Melville  et  al. 
(1997),  is  used  to  deform  the  grid  in  order  to  accommo¬ 
date  the  instantaneous  shape  of  the  membrane.  In  both 
situations,  great  care  must  be  exercised  to  avoid  grid- 
motion-induced  errors.  In  order  to  ensure  that  the  geo¬ 
metric  conservation  law  (GCL)  is  satisfied,  the  time 
metric  terms  are  evaluated  as  described  in  detail  in  Visbal 
and  Gaitonde  (2002). 

The  membrane  structural  equation  is  solved  numerically 
employing  the  procedure  discussed  in  Gordnier  (2008). 
Pinned  conditions  are  prescribed  at  the  membrane  leading- 
and  trailing-edges.  Coupling  of  the  aerodynamics  with  the 
membrane  response  occurs  through  the  imposed  pressure 
force,  A p,  in  Eq.  2  and  by  the  resulting  deflection  of  the 
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membrane,  which  is  returned  to  the  aerodynamic  grid.  A 
complete  synchronization  of  the  aerodynamic  and  struc¬ 
tural  systems  is  achieved  through  an  implicit  global 
subiteration  strategy  (Gordnier  2008).  During  each  subit¬ 
eration,  the  aerodynamic  loading  in  the  membrane 
equations  is  updated,  and  the  new  surface  displacements 
are  then  provided  to  the  aerodynamic  solver.  Using  this 
approach,  the  temporal  lag  between  the  aerodynamic  and 
structural  solutions,  as  well  as  numerical  linearization 
errors  may  be  effectively  eliminated. 

3  Transitional  flow  over  stationary  SD7003  airfoil 

In  this  section,  we  consider  the  effect  of  transition  on  the 
flow  structure  over  an  airfoil  at  Reynolds  numbers  relevant 
to  MAV  applications.  At  low  Reynolds  numbers 
( Rec  =  pcU^/p  <  105),  the  flow  may  remain  laminar  over 
a  significant  portion  of  the  airfoil,  and  is  unable  to  sustain 
even  mild  adverse  pressure  gradients.  For  moderate  inci¬ 
dence,  separation  leads  to  the  formation  of  a  laminar 
separation  bubble  (LSB)  which  breakdowns  into  turbulence 
prior  to  reattachment.  The  dynamics  of  this  laminar  sepa¬ 
ration  bubble,  as  it  moves  towards  the  leading  edge  with 
increasing  angle  of  attack  and  eventually  bursts,  has  a 
significant  impact  on  airfoil  performance. 

Results  from  large-eddy  simulations  for  an  SD7003 
airfoil  are  presented.  This  airfoil  geometry  was  selected 
since  it  exhibits  a  relatively  large  LSB  on  the  suction  side 
of  the  airfoil,  and  since  it  has  been  the  subject  of  recent 
experimental  and  computational  investigations  (Shyy  2008; 
01  et  al.  2005;  Radespiel  et  al.  2006;  Yuan  et  al.  2005; 
Galbraith  and  Visbal  2008).  Details  of  the  present  com¬ 
putations  can  be  found  in  Galbraith  and  Visbal  (2008)  and 
Galbraith  (2009)  wherein  the  effects  of  grid  resolution, 
Mach  number  and  computational  spanwise  extent  have 
been  explored.  Computations  were  performed  for  a  Rey¬ 
nolds  number  range  104  <  Rec  <9  x  104,  and  for  angles 
of  attack  2°  <  a  <  14°. 

3.1  Effect  of  angle  of  attack 

The  evolution  of  the  time-averaged  flow  structure  with 
increasing  incidence  is  examined  at  a  fixed  Reynolds 
number  ( Rec  =  6  x  104).  Shown  in  Fig.  1  are  the  mean 
streamline  patterns  and  the  Reynolds  stress  ( u'v ')  contours 
for  several  angles  of  attack.  The  corresponding  mean  sur¬ 
face  pressure  distributions  are  given  in  Fig.  2.  In  addition, 
Table  1  provides  a  summary  of  the  locations  of  laminar 
separation  and  turbulent  reattachment,  as  well  as  the 
maximum  bubble  height  (defined  in  terms  of  the  distance 
from  the  airfoil  surface  to  the  maximum  velocity  location 
at  the  edge  of  the  shear  layer).  Further  information  such  as 


skin-friction  coefficient  distribution  and  aerodynamic  loads 
can  be  found  in  Galbraith  and  Visbal  (2008). 

At  a  =  2°,  a  long  region  of  laminar  separation  is 
observed  over  the  airfoil  which  closes  just  upstream  of  the 
trailing  edge  ( xlc  =  0.92).  Relatively  small  values  of 
Reynolds  stress  magnitude  are  seen  near  the  reattachment 
zone.  As  the  angle  of  attack  increases  to  a  =  4°,  the  lam¬ 
inar  separation  point  moves  upstream  due  to  the  higher 
adverse  pressure  gradient  it  encounters.  The  reattachment 
location  is  also  displaced  upstream  due  to  the  more  intense 
process  of  transition  to  turbulence,  as  reflected  in  the 
increased  magnitude  of  Reynolds  stress.  As  a  result,  the 
mean  separation  bubble  contracts  both  in  height  and  in 
streamwise  extent.  This  overall  behavior  of  the  LSB  con¬ 
tinues  as  the  angle  of  attack  increases  up  to  a  =  11°.  Over 
this  range  of  incidence,  the  Cp  (Fig.  2)  distributions  display 
several  well-defined  characteristics,  including  a  suction 
peak,  a  plateau  associated  with  laminar  separation,  and  a 
fairly  rapid  recovery  following  transition.  All  of  these 
features  are  found  to  become  sharper  with  increasing  airfoil 
incidence.  At  a  =  14°,  the  laminar  separation  bubble  fails 
to  close  resulting  in  so-called  bubble  bursting.  Airfoil  stall 
is  evident  in  the  collapse  of  the  suction  peak  and  in  an 
essentially  flat  pressure  distribution  over  the  entire  suction 
surface.  A  large  mean  region  of  recirculation  (which 
extends  far  away  from  the  airfoil),  as  well  as  a  reduction  in 
the  maximum  magnitude  of  the  Reynolds  stress  are  also 
observed  (Fig.  1). 

Comparison  of  the  computed  time-averaged  flow  struc¬ 
ture  with  high-resolution  experimental  measurements  is 
shown  in  Fig.  3  for  a  =  4°.  The  measurements  denoted  as 
TU-BS  correspond  to  the  experiments  of  Radespiel  et  al. 
(2006),  whereas  the  data  referred  to  as  HFWT  was  obtained 
by  01  et  al.  (2005).  A  comparison  with  experiments  in 
terms  of  the  mean  separation  bubble  characteristics  is  also 
provided  in  Table  2.  Good  overall  agreement  is  observed 
between  the  experimental  and  computational  flow  struc¬ 
tures.  Some  discrepancies  in  specific  details  exist  between 
the  computations  and  experiment,  as  well  as  between  the 
experimental  data  sets  themselves.  In  particular,  the  LSB 
and  the  region  of  significant  Reynolds  stress  appear  to  be 
larger  in  the  computation  relative  to  experiment.  Also,  the 
reported  separation  locations  differ  between  the  two 
experiments,  due  in  part  to  the  difficulties  associated  with 
obtaining  the  precise  separation  point  from  near-wall  PIV 
measurements  for  a  very  shallow  separation  region.  It 
should  be  noted,  that  the  LSB  transition  process  is  expected 
to  be  quite  sensitive  to  several  factors  including  freestream 
turbulence  and  precise  experimental  setup  (e.g.,  aspect 
ratio,  interference  effects).  In  the  present  calculations,  no 
attempt  has  been  made  to  duplicate  any  such  factors  which 
in  general  are  not  well  characterized.  For  instance,  no 
freestream  turbulence  is  imposed  on  the  incoming 
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Fig.  1  Mean  streamline  pattern 
and  Reynolds  stress  as  a 
function  of  angle  of  attack  for 
Rec  =  6  x  104 
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computed  flow  fields,  which  may  account  for  differences  in 
the  LSB  size.  Elucidation  of  these  effects  would  require 
additional  computational  and  experimental  studies  wherein 
these  factors  could  be  varied  in  a  systematic  fashion. 
Nonetheless,  given  all  the  sources  of  uncertainty,  the  level 
of  agreement  observed  is  found  to  be  encouraging.  Addi¬ 
tional  comparison  with  experiments,  including  different 


angles  of  attack  and  aerodynamic  loads  is  provided  in 
Galbraith  and  Visbal  (2008). 

The  present  high-fidelity  simulations  permit  a  character¬ 
ization  of  the  complex  three-dimensional  flow  structure 
which  is  more  difficult  to  obtain  from  experimental  mea¬ 
surements.  The  evolution  of  the  computed  three-dimensional 
instantaneous  flow  structure  with  increasing  angle  of  attack 
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Fig.  2  Effect  of  angle  of  attack  on  mean  surface  pressure  distribu¬ 
tion,  Rec  =  6  x  104 


Table  1  Effect  of  angle  of  attack  on  computed  LSB  properties 


a  (degrees) 

Separation 

(xjc) 

Reattachment 

(x/c) 

Max  bubble 
height  (hb/c) 

2 

0.45 

0.92 

0.036 

4 

0.23 

0.65 

0.030 

6 

0.11 

0.45 

0.028 

8 

0.04 

0.28 

0.027 

11 

0.007 

0.16 

0.025 

14 

0.01 

- 

- 

Rec  =  6  x  104 


is  shown  in  Fig.  4.  This  figure  displays  iso-surfaces  of  the 
2-criterion  (Dubeif  and  Delcayre  2000)  which  provides  a 
means  of  visualizing  vortical  structures.  At  the  lowest 
incidence  (a  =  2°),  the  breakdown  of  the  separated  shear 
layer  begins  to  emerge  just  upstream  of  the  trailing  edge. 
Significant  changes  occur  with  a  small  increase  in  angle  of 
attack  to  a  =  4°.  Coherent  vortices  form  in  the  separated 
shear  layer  and  rapidly  breakdown  due  to  spanwise  insta¬ 
bility  effects.  Nonetheless,  the  original  coherent  spanwise 
structures,  surrounded  by  fine-scale  features,  are  still  dis- 
cernable.  With  further  increase  in  angle  of  attack  (up  to 
a  =  11°),  the  shear-layer  vortices  form  closer  to  the  leading 
edge,  and  the  viscous  flow  region  above  the  airfoil  is 
characterized  by  progressively  less  coherent  fine-scale 
features  associated  with  the  earlier  transition  to  turbulence. 
Finally,  for  a  =  14°,  the  shear  layer  fails  to  reattach  fol¬ 
lowing  separation  from  the  leading  edge  and  large  scale 
vortex  shedding  ensues.  These  results  demonstrate  that  the 


ILES  approach  is  capable  of  capturing  without  change  in 
parameters  the  entire  process  of  laminar  boundary  layer 
separation,  shear-layer  transition,  reattachment  of  the  LSB, 
and  eventually  the  passage  into  full  airfoil  stall. 

The  present  high-fidelity  3-D  simulations  indicate  that 
the  dynamics  of  the  flow  remains  qualitatively  unchanged 
prior  to  stall  (2°  <  a  <  11°).  In  all  cases,  the  flow  stays 
fairly  aligned  with  the  airfoil  surface,  and  transition  takes 
place  through  the  formation  of  a  closed  separation  bubble 
in  a  time-averaged  interpretation.  In  an  unsteady  sense, 
coherent  vortices  form  in  the  separated  shear  layer  and 
exhibit  spanwise  instabilities  and  subsequent  breakdown 
into  fine- scale  turbulent  structures.  High  magnitudes  of 
Reynolds  stresses  are  generated  which  promote  reattach¬ 
ment  to  the  surface.  This  process  becomes  more  abrupt  as 
the  angle  of  attack  and  the  corresponding  adverse  pressure 
gradient  increases,  resulting  in  a  shorter  mean  LSB.  This 
consistent  flow  evolution  for  a  stationary  airfoil  prior  to 
stall  therefore  appears  amenable  to  being  captured  by  2-D 
Reynolds-averaged  simulations  employing  simplified 
transition  models  (Shyy  2008;  Radespiel  et  al.  2006;  Yuan 
et  al.  2005).  Significant  challenges  however  are  anticipated 
when  using  those  approaches  for  stalled  flows  or  for  highly 
maneuvering  and  flexible  airfoils  involving  the  generation 
of  dynamic-stall-like  vortices. 

3.2  Effect  of  Reynolds  number 

In  order  to  explore  the  effect  of  Reynolds  number  on  the 
transition  process,  several  computations  were  performed 
for  a  fixed  angle  of  attack  (a  =  8°)  over  the  Reynolds 
number  range  104  <  Rec  <9  x  104.  Results  for  the  time- 
averaged  flow  fields  are  shown  in  Figs.  5,  6,  as  well  as  in 
Table  3. 

For  the  lowest  Reynolds  number  considered 
( Rec  =  104),  the  time-averaged  flow  exhibits  a  large 
recirculation  zone  above  the  airfoil.  The  corresponding 
surface  pressure  distribution  displays  a  very  weak  suction 
peak  and  a  fairly  flat  distribution  over  the  suction  side. 
Only  small  values  of  Reynolds  stress  are  observed  near  the 
trailing-edge  region.  For  Rec  =  3  x  104,  a  better  defined 
LSB  begins  to  emerge,  and  there  is  a  significant  increase  in 
the  magnitude  of  the  Reynolds  stress.  However,  the  LSB  is 
still  fairly  large  and  produces  significant  boundary-layer 
displacement.  The  corresponding  Cp  distribution  shows  a 
stronger  suction  peak  followed  by  a  plateau  and  a  mild 
recovery.  As  Reynolds  number  increases  further,  the  sep¬ 
aration  location  remains  effectively  unchanged  while  the 
LSB  shrinks  considerably  in  size.  This  is  also  accompanied 
by  a  more  effective  recovery  in  the  surface  pressure  dis¬ 
tribution.  The  variation  of  the  mean  drag  with  Reynolds 
number  for  a  fixed  angle  of  attack  is  given  in  Table  3.  It  is 
apparent  that  with  decreasing  Rec  there  is  a  significant 
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Fig.  3  Comparison  of 
computed  and  experimental 
mean  streamline  pattern, 
Reynolds  stress  and  velocity 
profiles  for  Rec  =  6  x  104 
and  a  =  4° 
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Table  2  Comparison  of  experimental  and  computed  LSB  properties 


Data  set 

Freestream 
turbulence  (%) 

Separation 
( xjc ) 

Reattachment 

(xjc) 

Max  bubble 
height  (hb/c) 

TU-BS 

0.08 

0.30 

0.62 

0.028 

HFWT 

0.1 

0.18 

0.58 

0.029 

ILES 

0 

0.23 

0.65 

0.030 

a  =  4°,  Rec  =  6  x  104 


increase  in  the  mean  drag  coefficient.  This  reduction  in 
aerodynamic  performance  is  a  direct  consequence  of  the 
significant  vertical  displacement  effect  created  by  the  lar¬ 
ger  LSB  at  lower  Reynolds  number  (Fig.  5). 

The  evolution  of  the  instantaneous  three-dimensional 
flow  structure  is  shown  in  Fig.  7,  using  again  an  iso-surface 
of  the  2-criterion.  For  Rec  >3  x  104,  the  process  of 
boundary  layer  transition  is  clearly  observed.  In  addition, 
finer  scale  turbulent  flow  features  are  generated  with 
increasing  Reynolds  number.  At  the  lowest  Reynolds 
number  considered  ( Rec  =  104),  shear  layer  transition  does 


not  occur  above  the  airfoil.  Instead,  the  flow  is  character¬ 
ized  by  the  formation  of  large-scale  coherent  vortices 
which  are  shed  from  the  leading  edge  and  which  only 
exhibit  mild  spanwise  instabilities  in  the  aft-region  of  the 
airfoil.  It  is  expected  that  the  effect  of  Reynolds  number  on 
the  transition  process  depends  also  on  angle  of  attack. 
Additional  results  showing  Reynolds  number  effects  for 
a  =  4°  and  a  =  14°  are  presented  in  Galbraith  (2009)  and 
Visbal  (2009).  At  the  lower  incidence  the  flow  failed  to 
transition  for  Rec  =  104.  In  the  stalled  regime  (a  =  14°), 
transition  of  the  separated  shear  layer  takes  place  in  the  aft- 
portion  of  the  airfoil  for  Rec  =  104. 

4  Transitional  flow  over  plunging  SD7003  airfoil 

In  this  section,  the  influence  of  a  prescribed  motion  on  the 
transitional  low-Reynolds-number  airfoil  flow  is  examined. 
Computations  were  performed  for  a  plunging  SD7003  air¬ 
foil  section  for  which  recent  PIV  measurements  (McGowan 
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Fig.  4  Evolution  of  instantaneous  three-dimensional  flow  structure  with  angle  of  attack,  Rec  =  6  x  104 


et  al.  2008)  have  been  conducted.  The  airfoil  is  set  at  a  static 
angle  of  attack  aG  =  4°,  and  the  reduced  frequency 
and  plunging  amplitude  are  k  =  nfc/U^  =  3.93  and 
h/c  =  0.05,  respectively.  These  motion  parameters  result  in 
a  maximum  excursion  in  induced  angle  of  attack  of  21.5° 
which  should  promote  leading-edge  separation  and  the 
generation  of  dynamic-stall-like  vortices.  Two  different 
Reynolds  numbers  ( Rec  =  104  and  4  x  104)  are  considered 
in  order  to  explore  the  impact  of  transition  on  the  flow 
structure  over  the  plunging  airfoil. 

Details  of  the  computational  mesh  generated  around  the 
SD7003  airfoil  are  described  elsewhere  (Visbal  2009).  Grid 
points  were  concentrated  near  the  airfoil  in  order  to  capture 
the  transition  process.  Therefore,  description  of  the  wake  is 
limited  to  the  region  just  downstream  of  the  airfoil  trailing 
edge  (up  to  x/c  =  2.0).  Study  of  the  downstream  evolution 
of  the  wake  far  away  from  the  airfoil  is  not  considered 
since  a  larger  grid  would  be  needed.  Both  2-D  and  3-D 
simulations  were  performed  in  order  to  explore  the  suit¬ 
ability  of  the  much  more  efficient  2-D  approach  for  low- 
Reynolds-number  applications,  as  well  as  to  compare  with 
previous  computations  (McGowan  et  al.  2008).  For  the  3-D 
calculations,  the  flows  were  assumed  to  be  periodic  in  the 


spanwise  direction  with  an  extent  sic  =  0.2.  This  value  of  s 
was  selected  based  on  a  previous  investigation  of  the  effect 
of  spanwise  extent  on  the  transitional  flow  past  the  same 
airfoil  geometry  at  fixed  angle  of  attack  (Galbraith  and 
Visbal  2008). 

Comparison  of  the  instantaneous  spanwise  vorticity 
component  obtained  from  both  2-D  and  3-D  computations 
is  shown  at  a  given  phase  of  the  plunging  motion  in  Fig.  8. 
For  Rec  =  104,  the  2-D  and  3-D  results  are  in  close 
agreement  with  each  other  over  a  significant  portion  of  the 
airfoil.  Some  discrepancies  exist  in  the  aft-portion  near  the 
trailing  edge,  as  well  as  in  the  near  wake  due  to  incipient 
3-D  effects.  However,  the  leading-edge  vortex  formation 
in  particular  exhibits  a  two  dimensional  character.  This 
behavior  appears  to  be  consistent  with  the  previous  sta¬ 
tionary  airfoil  (Sect.  3)  for  which  transition  was  not 
observed  at  this  low  Reynolds  number.  As  the  Reynolds 
number  is  increased  to  4  x  104,  significant  differences 
emerge  between  2-D  and  3-D  results.  The  coherent  vortices 
observed  in  the  2-D  simulations  breakdown  in  the  spanwise 
direction,  as  described  in  more  detail  later.  Given  the 
importance  of  spanwise  effects,  only  results  from  3-D 
computations  are  discussed  below. 
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Fig.  5  Effect  of  Reynolds 
number  on  time- averaged  flow 
structure,  a  =8° 
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■0.054 
-0  060 


Re60K 


uV/U3 

-0.005 
-0.011 
-0.017 
■0  023 
0  029 
-0.036 
■0.042 
-0  048 
-0.054 
-0  060 


uV/U* 

-Q  .005 
-0G11 
-0  017 
-0  023 
■0.029 
*0  QG8 
*0  042 
-0  048 
■0  054 
*0.060 


The  computed  instantaneous  flow  structure  over  the 
plunging  airfoil  for  Rec  =  104  is  shown  in  Fig.  9.  The 
selected  motion  phases  correspond  to  the  positions  of 
maximum  upward  displacement  (0  =  0),  maximum 
downward  velocity  (0  =  1/4),  maximum  downward  dis¬ 
placement  (O  =  1/2)  and  maximum  upward  velocity 
(0  =  3/4).  Although  not  included  here,  at  this  lower 
Reynolds  number,  the  overall  instantaneous  and  phase- 
averaged  flow  structures  are  found  to  be  in  reasonable 
agreement  with  each  other  due  to  the  limited  impact  of 
transition  (Visbal  2009).  A  fairly  good  comparison 
between  experimental  and  phase-averaged  vorticity  is  also 
shown  in  Visbal  (2009). 


Due  to  the  large  angle  attack  induced  by  the  plunging 
motion  during  the  downstroke,  leading-edge  separation 
occurs  on  the  upper  surface  of  the  airfoil,  and  coherent 
dynamic- stall-like  vortices  are  generated.  Two  distinct 
leading-edge  vortices  (denoted  as  ‘1’  and  ‘2’  in  Fig.  9c)  are 
formed  by  the  time  the  airfoil  reaches  its  bottom  dead 
center.  During  the  upstroke,  these  vortices  propagate  close 
to  the  airfoil  surface  thereby  precluding  massive  stall  for 
this  high  value  of  reduced  frequency.  Given  the  short 
period  of  the  motion  (T  =  n/k  =  0.8),  the  dynamic  stall 
vortices  emerge  before  the  pair  of  vortices  generated  in  the 
previous  cycle  (denoted  as  ‘3’  and  ‘4’  in  Fig.  9d)  have 
reached  the  airfoil  trailing  edge.  In  addition  to  the  primary 
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Fig.  6  Effect  of  Reynolds  number  on  mean  surface  pressure  distri¬ 
bution,  a  =  8° 


Table  3  Effect  of  Reynolds  number  on  computed  LSB  propreties 
and  mean  drag  coefficient 


Rec 

Separation 
( xjc ) 

Reattachment 

(Xr/C) 

Max  bubble 
height  (hb/c) 

Mean  drag 
(CflmeJ 

1  X 

104 

0.09 

0.98 

0.217 

0.082 

3  x 

104 

0.05 

0.53 

0.073 

0.070 

6  x 

104 

0.04 

0.28 

0.027 

0.043 

9  x 

104 

0.04 

0.20 

0.014 

0.035 

a  =  8° 


leading-edge  vortical  structures,  another  distinct  feature  is 
the  ejection  of  vorticity  of  the  opposite  sign  due  to  the 
ensuing  vortex  surface  interaction.  This  ejected  vorticity  is 
quite  prominent  between  the  two  primary  dynamic-stall 
vortices  (Fig.  9d),  and  eventually  surrounds  completely  the 
leading  vortex  (denoted  as  ‘1’).  Further  downstream,  this 
secondary  vorticity  becomes  less  apparent  due  to  spanwise 
instability  effects,  as  discussed  below.  Separation  and 
formation  of  a  single  dynamic-stall  vortex  is  also  observed 
on  the  airfoil  lower  surface  (Fig.  9b)  as  a  result  of  the  large 
negative  angle  of  attack  induced  during  the  upstroke. 

The  present  implicit  large-eddy  simulations  permit  a 
detailed  description  of  the  3-D  instantaneous  flow  structure 
which  complements  experimental  planar  observations.  The 
overall  instantaneous  3-D  flow  features  for  Rec  =  104  are 
displayed  for  several  phases  in  Fig.  10.  These  plots  show 
an  iso- surface  of  vorticity  magnitude  colored  by  density  in 
order  to  enhance  contrast.  At  all  phases  of  the  plunging 
motion  shown,  one  can  observe  coherent  spanwise  vortical 


structures  or  rollers  which  exhibit  a  fairly  2-D  character. 
These  vortices  are  surrounded  by  complex  3-D  flow  fea¬ 
tures  resulting  from  spanwise  instabilities.  In  particular,  the 
dynamic  stall  vortex  forming  near  the  leading  edge  retains 
its  spanwise  coherence;  however,  the  vorticity  ejected  from 
the  airfoil  surface  (as  a  result  of  the  strong  vortex-surface 
interaction)  rapidly  breakdowns  giving  rise  to  the  appear¬ 
ance  of  longitudinal  vortical  structures.  Therefore,  for  this 
high  reduced  frequency  and  low  Reynolds  number,  tran¬ 
sitional  effects  over  the  airfoil  appear  to  be  minor,  and  are 
associated  with  the  ejected  secondary  vorticity.  However, 
as  discussed  next,  this  situation  changes  significantly  with 
increasing  Reynolds  number. 

Results  for  Rec  =  4  x  104  are  considered  next.  Fig¬ 
ure  11  displays  a  comparison  of  the  experimental  and 
computed  phase-averaged  streamwise  velocity  and  span- 
wise  vorticity  contours  at  the  end  of  the  upstroke.  The 
computed  results  are  found  to  be  in  reasonable  qualitative 
agreement  with  the  experimental  measurements.  A  more 
quantitative  comparison  is  shown  in  Fig.  12  in  terms  of  the 
streamwise  velocity  profiles  in  the  near  wake  at  xlc  =1.5. 
The  agreement  between  the  measured  and  computed 
profiles  is  in  general  quite  good.  The  velocity  profile  cor¬ 
responding  to  the  phase  of  maximum  vertical  displacement 
(0  =  0)  exhibits  a  well-defined  jet  character.  This  jet-like 
behavior  is  consistent  with  the  fact  that  the  trailing-edge 
vortex  shed  from  the  airfoil  lower  surface  is  above  the 
vortex  shed  from  the  upper  surface  boundary  layer  (see 
Fig.  11c,  d).  As  it  will  be  discussed  later,  this  inverted 
vortex  street  results  in  a  net  mean  thrust. 

The  computed  instantaneous  and  phased-averaged 
spanwise  vorticity  fields  are  shown  in  Fig.  13.  Significant 
differences  are  observed  between  the  instantaneous  and 
phased-averaged  representations  which  reflect  the  fully  3-D 
transitional  nature  of  the  flow.  The  instantaneous  vorticity 
exhibits  fine-scale  flow  features  within  the  dynamic  stall 
vortices,  as  well  as  in  the  near  wake.  These  features  are 
essentially  eliminated  in  the  phase-averaging  process. 
Increasing  the  Reynolds  number  has  clearly  modified  the 
details  of  the  phase-averaged  vortical  structure  emerging 
from  the  unsteady  separation  process  near  the  airfoil 
leading  edge.  As  discussed  earlier,  for  Rec  =  104  (Fig.  9c, 
d),  the  leading-edge  vortex  system  exhibited  a  pair  of 
dynamic- stall-like  vortices  which  retained  their  identity  as 
they  propagated  downstream  in  the  absence  of  significant 
spanwise  transitional  effects.  By  contrast,  for  Rec  =  4  x 
104  (Fig.  13),  a  single  and  more  diffused  leading-edge 
vortex  is  observed  in  the  phased-averaged  representation. 

The  computed  instantaneous  3-D  flow  structure  above 
the  plunging  airfoil  for  Rec  =  4  x  104  is  displayed  in 
Fig.  14.  The  most  prominent  overall  flow  features  are  the 
leading-edge  vortices  and  their  breakdown  into  fine-scale 
turbulence  due  to  spanwise  instability  effects.  At  the  end  of 
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Fig.  7  Effect  of  Reynolds  number  on  three-dimensional  instantaneous  flow  structure,  a  =  8° 


Fig.  8  Comparison  of 
instantaneous  spanwise  vorticity 
for  2-D  and  3-D  plunging  airfoil 
simulations 


the  upstroke  (O  =  0),  the  boundary  layer  near  the  leading 
edge  remains  laminar  and  attached  to  the  airfoil  surface.  By 
the  time  the  airfoil  has  reached  its  maximum  downward 
velocity  (O  =  1/4),  a  small  region  of  laminar  separation 
has  emerged  near  the  leading  edge.  Between  0=1/4  and 
O  =  1/2,  this  separation  region  evolves  into  a  coherent 
dynamic- stall-like  vortex  which  subsequently  breakdowns 
into  fine-scale  structures.  This  transitional  leading-edge 


vortex  propagates  downstream  close  to  the  airfoil  surface, 
as  seen  in  Fig.  13.  Given  the  relatively  high  value  of 
reduced  frequency,  the  dynamic- stall-like  vortex  generated 
in  the  previous  cycle  (denoted  as  V\  in  Fig.  14c)  has  not 
reached  the  airfoil  trailing  edge  before  a  new  leading-edge 
vortex  (V2)  is  generated. 

The  transition  of  the  dynamic  stall  vortex  near  the 
leading  edge  is  described  in  detail  in  Visbal  (2009).  This 
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Fig.  9  Instantaneous  spanwise 
vorticity  at  selected  phases  of 
the  plunging  motion  for 
Rec  =  104 


Fig.  10  Iso- surface  of 
instantaneous  vorticity 
magnitude  at  selected  phases 
of  the  plunging  motion  for 
Rec  =  104 


Fig.  11  Comparison  of 
computed  and  experimental 
phased-averaged  streamwise 
velocity  (a,  b)  and  spanwise 
vorticity  (c,  d)  at  the  end  of  the 
upstroke  for  Rec  =  4  x  104 


<C) 


Experiment 

(b) 


0.0  0,2  0.4  0.6  0.3  1.0  12  14 


(d) 


transition  process  is  found  to  be  rather  abrupt  and  takes 
place  over  a  short  non-dimensional  time  interval  of  order 
AtU^/c  <0.04.  This  represents  a  fast  onset  compared  to 
either  the  plunging  motion  or  the  vortex  convection  time 
scales.  If  remains  to  be  seen  whether  this  process  can  be 
duplicated  by  standard  Reynolds-averaged  eddy-viscosity 
models  which  rely  on  dissipation-like  terms.  In  addition, 


the  mixed  laminar-transitional  character  of  this  highly 
unsteady  flow  is  quite  apparent.  For  instance,  the  vortex 
structure  of  Fig.  14c  displays  an  extensive  region  of  fairly 
laminar  flow  between  the  two  transitional  leading-edge 
vortices  generated  in  consecutive  cycles  of  the  plunging 
motion.  This  mixed  flow  structure,  although  amenable  to 
simulation  with  the  present  high-fidelity  ILES  approach, 
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Fig.  12  Comparison  of 
computed  and  experimental 
phased-averaged  streamwise 
velocity  profiles  in  the  near 
wake  (x/c  =  1.5)  of  the 
plunging  airfoil  for 
Rec  =  4  x  104 


(a)  ®  =  o 


(b)  ©■  1/4 


Fig.  13  Comparison  of 
computed  instantaneous  and 
phased-averaged  spanwise 
vorticity  at  selected  phases  of 
the  plunging  motion  for 
Rec  =  4  x  104 
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will  remain  a  significant  challenge  for  more  efficient 
Reynolds-averaged  methods  incorporating  simplified  tran¬ 
sition  models. 

The  time  histories  of  the  lift  coefficient  for  both  Rey¬ 
nolds  numbers  are  shown  in  Fig.  15a.  Clearly,  CL  appears 
to  be  essentially  independent  of  Reynolds  number  and  of 
the  3-D  transitional  aspects  of  the  flow  field.  The  insensi¬ 
tivity  of  the  aerodynamic  lift  is  due  to  the  fact  that  at  this 
high  value  of  reduced  frequency,  CL  is  dominated  by  the 
acceleration  of  the  airfoil  which  scales  with  k 2 .  Indeed,  the 


computed  lift  is  found  to  be  in  close  agreement  with  the 
prediction  given  by  the  in  viscid  theory  (Fung  1993) 
(Fig.  15b)  which  in  turn  is  shown  to  be  dominated  by  non- 
circulatory  effects.  The  computed  mean  drag  coefficients 
for  Rec  =  104  and  Rec  =  4  x  104  were  found  to  be 
—0.055  and  —0.083,  respectively,  where  negative  CD  val¬ 
ues  correspond  to  net  thrust  on  the  airfoil.  As  it  is  perhaps 
to  be  expected,  the  drag  coefficient  shows  a  significant 
dependence  on  Reynolds  number  and  on  the  transitional 
aspects  of  the  flow. 
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Fig.  14  Iso- surface  of 
instantaneous  vorticity 
magnitude  showing  flow 
structure  above  plunging  airfoil 
at  selected  phases  of  the 
plunging  motion  for 
Rec  =  4  x  104 


Fig.  15  Computed  lift 
coefficient  histories  for 
plunging  airfoil 


5  Flexible  membrane  airfoil 

In  this  section,  the  impact  of  fluid-structure  interaction  on 
the  low-Reynolds-number  flow  past  a  flexible  airfoil  is 
investigated.  The  membrane  airfoil  configuration  consid¬ 
ered  is  based  on  the  experimental  model  of  Rojratsirikul, 
Wang,  and  Gursul  (Rojratsirikul  et  al.  2008).  In  these 
experiments,  the  membrane  wing  consists  of  a  thin  latex 
sheet  stretched  between  two  small,  aerodynamically- 
shaped,  rigid  mounts.  The  sheet  is  glued  to  the  mounts 
resulting  in  pinned  boundary  conditions  for  the  flexible 
membrane.  Experimental  observations  (Rojratsirikul  et  al. 
2008)  indicate  that  the  membrane  deformation  is  essen¬ 
tially  uniform  in  the  spanwise  direction.  This  observation 
justifies  the  initial  assumption  of  2-D  deformation,  at  least 
for  low  Reynolds  number  simulations. 

For  the  present  computations,  a  freestream  Mach  num¬ 
ber  Mqo  =  0-05  is  prescribed  in  order  to  approach  the 
incompressible  situation.  Unless  otherwise  noted,  the 
results  presented  are  for  a  Reynolds  number  based  on  air¬ 
foil  chord  of  Rec  =  2,500.  The  static  angle  of  attack  ranged 


from  a  =  4°  to  a  =  20°.  The  membrane  structural 
parameters  in  Eq.  2  are  specified  as  follows:  mass  ratio, 
psh  =  0.589,  membrane  rigidity,  Eh  =  50,  no  structural 
damping,  Cd  =  0,  and  no  membrane  pretension,  S0  =  0. 
To  explore  the  impact  of  transition  an  additional  compu¬ 
tation  at  a  higher  Reynolds  number  ( Rec  =  4.85  x  104), 
and  with  structural  parameters  Eh  =  14.48  and  psh  =1.2 
is  also  presented.  Further  simulations  for  higher  Reynolds 
numbers  and  for  other  values  of  the  structural  parameters 
are  described  in  Gordnier  (2008)  and  Gordnier  and  Attar 
(2009)  wherein  additional  details  of  the  computations  can 
also  be  found. 

The  evolution  of  the  time-averaged  flow  structure  with 
increasing  angle  of  attack  is  shown  in  Figs.  16,  17,  and  18 
in  terms  of  the  mean  vorticity  and  streamline  patterns,  the 
membrane  shape  and  pressure  distributions,  as  well  as  the 
mean  aerodynamic  lift  coefficient.  The  most  obvious  effect 
of  flexibility  is  the  development  of  significant  camber  in 
the  mean  shape  of  the  membrane  due  to  the  differences  in 
pressure  between  the  upper  and  lower  sides  of  the  airfoil. 
The  maximum  deflection  (Fig.  17a)  grows  over  the 
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Fig.  16  Mean  vorticity  and 
streamline  patterns  for 
membrane  airfoil  at  selected 
angles  of  attack:  a,  f  a  =  4°;  b, 
g  a  =  8°;  c,  h  a  =  12°;  (d,  i) 
a  =  16°;  and  (e,  j)  a  =  20° 


incidence  range  considered  from  zmax/c  =  0.0724  at 
a  =  4°  to  ZmaJc  =  0.0899  at  a  =  20°.  For  a  =  4°  the 
membrane  shape  is  nearly  symmetric.  However,  it  becomes 
progressively  more  asymmetric  with  increasing  incidence. 
For  a  =  16°,  the  point  of  maximum  deflection  reaches 
xlc  =  0.426  before  it  starts  moving  back  downstream  with 
further  increase  in  angle  of  attack. 

The  impact  of  increasing  angle  of  attack  (and  corre¬ 
spondingly  mean  camber)  on  the  overall  mean  flowfield 
is  shown  in  Fig.  16.  At  a  =  4°,  the  flow  detaches  at 
xlc  =  0.738  resulting  in  a  narrow  separation  bubble  which 
extends  to  the  trailing  edge  (Fig.  16f).  A  very  small  sepa¬ 
ration  bubble  is  present  near  the  leading  edge  associated 
with  the  rigid  membrane  mount.  It  should  be  noted  that 
separation  was  always  found  to  occur  at  the  rigid  mount. 
This  fact  suggests  that  how  the  underlying  support  structure 
for  a  membrane  is  designed  may  have  direct  influence  on 
the  flexible  airfoil  performance.  As  the  angle  of  attack  is 
raised  to  a  =  8°,  the  downstream  separation  point  moves 
forward  to  x  =  0.57,  and  the  normal  extent  of  the  separa¬ 
tion  bubble  increases.  The  leading-edge  separation  bubble 
also  increases  in  size  with  attachment  occurring  at 
x  =  0.161.  By  a  =  12°  the  two  separation  zones  have 
merged  creating  one  large  separation  region  that  is  narrow 
in  vertical  extent  and  covers  nearly  the  full  length  of  the 
airfoil.  At  a  =  16°,  the  separation  zone  has  shrunk  with 
attachment  at  x  =  0.774,  but  the  normal  extent  of  the 
separated  region  has  grown.  Finally,  by  a  =  20°,  a 


massively  separated  flow  region  exists  over  the  airfoil.  A 
small  secondary  separation  zone  is  also  observed  near  the 
leading  edge  of  the  membrane. 

The  corresponding  evolution  of  the  mean  vorticity  field 
is  shown  in  Fig.  16a-e.  At  the  lowest  angle  of  attack 
(Fig.  16a),  the  shear  layer  separates  from  the  upstream 
mount  but  rapidly  attaches  to  the  membrane  airfoil.  It 
subsequently  departs  from  the  membrane  surface  down¬ 
stream  at  the  previously  noted  separation  location.  By 
a  =  12°  (Fig.  16c),  the  shear  layer  separates  from  the 
upstream  mount  and  fails  to  reattach  to  the  membrane 
surface.  A  significant  thickening  of  the  shear  layer  is  also 
observed.  At  the  higher  angles  of  attack,  the  departure 
angle  of  the  shear  layer  from  the  leading  edge  mount  grows 
and  the  shear  layer  moves  increasingly  away  from  the 
membrane.  At  a  =  20°  (Fig.  16e),  a  small  region  of  vor¬ 
ticity  of  the  opposite  sign  is  located  under  the  shear  layer 
near  the  leading  edge  associated  with  the  development  of 
secondary  separation. 

The  resulting  mean  surface  pressure  distributions  on  the 
membrane  are  plotted  in  Fig.  17b.  As  the  incidence  angle 
increases,  a  well-defined  suction  peak  develops  on  the 
upper  surface  in  the  leading  edge  region.  On  the  lower 
surface,  the  pressure  progressively  increases  with  angle  of 
attack. 

Shown  in  Fig.  18  is  a  comparison  of  the  mean  lift 
coefficient  for  the  flexible  membrane  airfoil,  as  well  as  for 
a  rigid  flat  airfoil.  The  CL  distribution  for  a  symmetric 
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Fig.  17  Mean  membrane 
deflection  (a)  and  surface 
pressure  coefficient  (b)  for 
various  angles  of  attack 


Fig.  18  Comparison  of  mean  lift  coefficient  for  flexible  membrane 
and  for  rigid  symmetric  and  cambered  wing  configurations 


airfoil  based  on  thin-airfoil  theory  is  also  shown  for  ref¬ 
erence  purposes.  The  rigid  airfoil  lift  coefficient  increases 
approximately  in  a  linear  fashion  with  a  slope  which  is 
close  to  the  theoretical  value.  A  maximum  CL  is  achieved 
for  a  =  12°.  Beyond  this  angle  of  attack,  the  lift  distribu¬ 
tion  displays  the  characteristics  of  a  stalled  flow.  The 
flexible  membrane  airfoil  exhibits  a  lift  variation  which 
differs  significantly  from  its  rigid  counterpart.  If  one 
extrapolates  the  lift  coefficient  to  zero  incidence,  a  signif¬ 
icant  lift  is  achieved  due  to  the  camber  effect  induced  by 
flexibility.  When  the  angle  of  attack  is  increased  up  to 
a  =  12°,  CL  increases  in  an  approximate  linear  fashion 
albeit  with  a  reduced  slope  relative  to  either  the  rigid  case 
or  the  theoretical  distribution.  For  higher  incidence 
(12°  <  a  <  20°),  the  flexible  airfoil  lift  coefficient  contin¬ 
ues  to  rise,  although  there  is  a  further  reduction  in  the  lift- 
curve  slope.  Therefore,  at  high  angles  of  attack,  airfoil 


flexibility  results  in  substantial  lift  enhancement  relative  to 
the  rigid  case.  For  the  range  of  a  considered,  there  appears 
to  be  no  evidence  of  stall  for  the  membrane  airfoil,  at  least 
in  terms  of  the  aerodynamic  lift  coefficient.  Lift  enhance¬ 
ment  due  to  wing  flexibility  at  high  angles  of  attack  has 
also  been  demonstrated  in  wind  tunnel  tests  by  Waszak 
et  al.  (2001)  and  is  due  to  both  the  effect  of  the  induced 
mean  camber,  as  well  as  the  dynamic  surface  motion  to  be 
described  next. 

Instantaneous  snapshots  of  the  vorticity  field  (Fig.  19) 
illustrate  the  evolving  unsteady  flow  structure  with 
increasing  angle  of  attack.  At  a  =  4°,  the  flow  separates 
just  past  the  midpoint  of  the  airfoil.  The  separated  shear 
layer  does  not  become  unstable  or  roll  up  until  it  passes  the 
trailing  edge  and  interacts  with  the  separating  flow  from  the 
lower  surface.  Therefore,  the  flow  over  the  airfoil  is  pre¬ 
dominantly  steady.  As  the  angle  of  attack  is  increased, 
the  onset  of  the  shear  layer  rollup  moves  progressively 
upstream  and  the  flow  becomes  increasingly  unsteady  over 
the  airfoil.  Up  to  a  =  12°,  the  vortices  remain  close  to  the 
airfoil  surface  as  they  propagate  downstream.  At  a  =  16°, 
the  vortex  train  begins  to  depart  from  the  membrane.  The 
convecting  shear-layer  vortices  interact  with  the  surface 
boundary  layer  resulting  in  the  development  of  secondary 
vortical  structures,  and  the  eventual  eruption  of  vorticity  of 
the  opposite  sign  from  the  surface.  Pairing  of  the  shear- 
layer  vortices  is  also  observed  in  computed  animations  for 
the  higher  angles  of  incidence. 

The  unsteady  structural  response  of  the  membrane  air¬ 
foil  is  shown  in  Fig.  20.  This  figure  displays  the  membrane 
oscillations  by  means  of  x—t  diagrams  showing  the  dif¬ 
ference  between  the  instantaneous  and  mean  membrane 
deflections.  The  membrane  is  nominally  stationary  at 
a  =  4°  since,  as  noted  earlier,  the  flow  is  steady  over  the 
major  portion  of  the  airfoil  (Fig.  19a).  At  a  =  8.0°,  a 
third  mode  standing  wave  response  develops  with  a 
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(a) 


Fig.  19  Instantaneous  vorticity  for  membrane  airfoil  at  selected 
angles  of  attack:  a  a  =  4°,  b  a  =  8°,  c  a  =  12°,  d  a  =  16°,  and  e 
a  =  20° 

nondimensional  frequency,  St  =  fc/u^  =  1.43.  This  fre¬ 
quency  correlates  within  the  tolerance  of  the  Fourier 
analysis  to  the  vortex  shedding  frequency  St  «  1.45  as 
measured  at  a  point  above  the  airfoil  by  which  the  shed 
vortices  are  passing.  At  the  higher  angles  of  attack  a  =  12° 
to  a  =  20°  the  membrane  structural  response  exhibits  a 
less  regular  behavior  resulting  from  a  combination  of 
structural  modes.  Multiple  frequency  peaks  are  found  in  the 
Fourier  analysis  of  the  structural  response.  By  a  =  20° 
three  peaks  emerge  in  the  spectral  analysis  of  the  structural 
response  with  the  dominant  frequency  being  St  =  0.59. 
This  corresponds  to  the  computed  vortex  shedding  fre¬ 
quency,  St  =  0.6,  which  correlates  with  the  standard  value 
for  the  wake  shedding  frequency  when  rescaled  by  angle  of 
attack,  St  sin(a)  =  0.2.  The  maximum  peak-to-peak 
amplitude  of  the  deflection  reached  values  of  order  0.07  c 
which  is  of  the  same  order  as  the  maximum  mean  deflec¬ 
tion  (Fig.  17a).  Therefore,  in  addition  to  the  mean  camber 
effect  previously  discussed,  flexibility  provides  a  signifi¬ 
cant  forcing  of  the  unsteady  separated  flow  over  the 
membrane. 


In  order  to  separate  the  effects  of  the  dynamic  motion  of 
the  membrane  from  the  mean  induced  camber  at  higher 
angles  of  attack,  a  computation  was  run  at  a  =  20°  for  a 
rigid  cambered  airfoil  with  its  shape  prescribed  by  the 
previously  computed  mean  membrane  deflection.  A  com¬ 
parison  of  the  mean  and  unsteady  solutions  for  these  two 
cases  is  shown  in  Fig.  21.  The  mean  streamline  patterns 
(Figs.  21a,  b)  show  that  the  rigid  airfoil  flowfield  exhibits  a 
larger  recirculation  region,  as  well  as  a  larger  secondary 
separation  bubble  located  further  downstream  on  the  air¬ 
foil.  These  significant  differences  between  the  time- 
averaged  flows  for  the  rigid  and  flexible  airfoils  result  from 
a  distinct  change  in  the  instantaneous  vortical  structure 
(Fig.  21c,  d).  The  dynamic  motion  of  the  membrane  excites 
the  shear  layer  separating  at  the  leading  edge  causing  it  to 
roll  up  sooner  and  to  form  a  series  of  more  coherent  vor¬ 
tices  which  convect  closer  to  the  airfoil  surface.  The 
trailing  edge  vortex  in  the  dynamic  case  is  also  reduced  in 
strength  and  forms  farther  downstream  and  away  from  the 
airfoil.  A  similar  comparison  between  a  rigid  cambered 
airfoil  and  a  flexible  membrane  has  been  investigated 
experimentally  by  Rojratsirikul  et  al.  (2008).  Their  flow 
visualizations  (Rojratsirikul  et  al.  2008)  also  demonstrate  a 
reduction  of  the  separation  zone  due  to  the  dynamic  motion 
of  the  membrane.  The  overall  impact  of  flexibility,  through 
the  combined  effects  of  induced  mean  camber  and  dynamic 
motion,  is  a  significant  delay  in  the  stall  of  the  airfoil  with  a 
corresponding  increase  in  lift  coefficient  in  the  post-stall 
region  (Fig.  18).  This  enhanced  lift  arises  predominantly 
from  the  effect  of  mean  camber  with  additional  lift 
enhancement  resulting  from  the  membrane  motion.  Further 
investigation  of  the  impact  of  flexibility  on  the  overall 
aerodynamics  is  warranted  as  it  may  provide  a  means  for 
passively  controlling  the  flow  in  order  to  achieve  improved 
airfoil  performance  through  judicious  aeroelastic  tailoring. 

As  the  Reynolds  number  is  increased  from  the  low 
values  considered  to  this  point  to  more  moderate  values, 
the  impact  of  transition  on  the  flow  becomes  important. 
Figure  22  compares  a  laminar,  2-D  simulation  with  a  3-D 
implicit  LES  simulation  for  Rec  =  4.85  x  104  and  a  =  8°. 
For  the  2-D  laminar  simulation  where  transition  of  the  flow 
cannot  be  captured  directly  (Fig.  22a),  strong  vortices  are 
formed  and  shed  from  the  leading  edge.  As  these  vortices 
convect  downstream,  they  engender  a  very  dynamic 
interaction  with  the  flexible  membrane  resulting  in  large 
higher- mode  fluctuations.  The  results  of  the  3-D  compu¬ 
tation  using  the  implicit  LES  approach  (Fig.  22b),  exhibit 
rapid  transition  of  the  flow  separating  from  the  leading 
edge,  and  the  subsequent  development  of  large-scale 
coherent  turbulent  structures  with  embedded  fine  scale 
features.  A  more  regular  second-mode  oscillation  of  the 
membrane  results  for  this  transitional  flow.  It  is  therefore 
apparent  that  correct  prediction  of  the  transitional  behavior 
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Fig.  20  x—t  diagram  of 
instantaneous  membrane 
deflection  for  selected  angles  of 
attack:  a  a  =  4°,  b  a  =  8°,  c 
a  =  12°,  d  a  =  16°,  and  e 
a  =  20° 
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Fig.  21  Comparison  of  rigid 
cambered  airfoil  (a,  c)  and 
flexible  membrane  (b,  d)  flow 
fields  for  a  =  20°:  a,  b  mean 
streamlines,  and  c,  d 
instantaneous  spanwise  vorticity 


is  critical  in  order  to  capture  the  proper  membrane  struc¬ 
tural  response.  Therefore,  computational  approaches  that 
are  capable  of  treating  these  unsteady,  mixed  transitional 
flows  are  required  at  these  moderate  Reynolds  numbers. 
Further  discussions  of  the  ILES  simulations  for  membrane 
airfoils  may  be  found  in  Gordnier  and  Attar  (2009). 

6  Summary  and  conclusions 

A  high-fidelity  implicit  large-eddy  simulation  approach  has 
been  employed  to  investigate  several  aspects  of  the  com¬ 
plex  unsteady  flow  physics  encountered  by  small-scale 
vehicles  in  low-speed  flight.  This  regime  is  characterized 


by  highly  unsteady  mixed  laminar-transitional-turbulent 
flows  over  moving  and  flexible  surfaces.  Three  separate 
fluid  dynamic  phenomena  are  considered,  including:  lam¬ 
inar  separation  and  transition  over  a  stationary  airfoil, 
transition  effects  on  the  dynamic  stall  vortex  generated  by  a 
plunging  airfoil,  and  the  effect  of  flexibility  on  the  flow 
structure  above  a  membrane  airfoil.  Comparison  with  high- 
resolution  PIV  experiments  has  been  provided  whenever 
possible. 

The  process  of  boundary-layer  transition  over  a  sta¬ 
tionary  SD7003  airfoil  has  been  examined  for  a  range  of 
Reynolds  number  (104  <  Rec  <9  x  104)  and  angle  of 
attack  (2°  <  a  <  14°).  Prior  to  stall,  the  flow  is  character¬ 
ized  by  a  separated  shear  layer  which  exhibits  the 
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Fig.  22  Effect  of  transition  at 
Reynolds  number, 

Rec  =  48,500  :  a  2-D  laminar 
flow  simulation,  b  3-D 
transitional  flow  simulation 


formation  of  span  wise  coherent  vortical  structures.  These 
vortices  breakdown  due  to  spanwise  instability  mecha¬ 
nisms  and  promote  the  reattachment  of  the  boundary  layer. 
In  a  time-averaged  representation,  a  laminar  separation 
bubble  (LSB)  is  formed  which  diminishes  both  in  stream- 
wise  and  vertical  extent  with  increasing  angle  of  attack.  At 
fixed  incidence,  decreasing  Reynolds  number  promotes  the 
formation  of  a  larger  LSB.  By  Rec  =  104,  transition  fails  to 
occur  over  the  airfoil  with  an  accompanying  significant 
increase  in  drag  coefficient.  The  sensitivity  of  these  tran¬ 
sitional  flows  to  freestream  turbulence  or  other  imposed 
disturbances  needs  to  be  elucidated. 

The  unsteady  low-Reynolds-number  ( Rec  =  104  and 
4  x  104)  flow  past  a  SD7003  airfoil  plunging  at  a  relatively 
high  reduced  frequency  ( k  =  3.93)  was  investigated. 
Comparison  was  made  with  recent  PIV  measurements.  For 
the  parameters  chosen,  the  large  motion-induced  angle  of 
attack  generates  dynamic-stall-like  vortices  in  the  leading- 
edge  region  both  above  and  underneath  the  airfoil.  Due 
to  the  high  value  of  k  and  small  plunging  amplitude 
(h/c  =  0.05),  these  vortices  convect  downstream  close  to 
airfoil  surface  precluding  massive  stall.  For  Rec  =  104,  the 
computed  instantaneous  and  phased-averaged  flow  struc¬ 
tures  are  in  fairly  good  agreement  with  each  other.  This  is 
due  to  the  incipient  process  of  transition  which  is  mainly 
limited  to  the  secondary  vorticity  ejected  from  the  airfoil 
surface.  When  the  Reynolds  number  is  increased  to 
4  x  104,  the  dynamic-stall  vortices  breakdown  shortly 
after  their  onset  due  to  spanwise  instability  effects.  The 
instantaneous  flow  exhibits  fine- scale  structures  which  are 
eliminated  in  the  phase-averaged  representation.  At  this 
intermediate  Reynolds  number,  a  quite  complex  mixed 
transitional-turbulent  flow  field  exists  which  would  pose  a 
formidable  challenge  for  its  simulation  employing  simpli¬ 
fied  approaches. 


The  effect  of  structural  compliance  on  the  unsteady  flow 
past  a  membrane  airfoil  at  low  Reynolds  number 
(Rec  =  2.5  x  103)  has  been  explored  by  solving  the  cou¬ 
pled  fluid- structural  system.  The  most  notable  effect  of 
membrane  flexibility  is  the  introduction  of  a  mean  camber 
due  to  the  fluid  forces  acting  on  both  sides  of  the  airfoil.  At 
high  incidence,  large  membrane  fluctuations  are  also 
established  with  peak-to-peak  amplitudes  commensurate 
with  the  mean  deflection.  This  dynamic  behavior  provides 
an  aeroelastically  induced  excitation  of  the  separated  shear 
layer  which  rolls  up  more  readily  and  remains  closer  to  the 
membrane.  These  combined  effects  result  in  an  increase  in 
lift  and  in  a  delay  in  the  onset  of  stall.  At  a  higher  Reynolds 
number  (Rec  =  4.85  x  104),  the  inclusion  of  transition 
effects  is  required  to  capture  the  correct  membrane  struc¬ 
tural  response. 
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